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Abstract. The validity of the classical formula for the rate of secular apsidal motion in close binaries is investigated 
for a sequence of models of a 5 Mq star ranging from the last stages of the C^'^ — > N'^* reaction to the phase where 
hydrogen is exhausted in the core. For binaries with short orbital periods, the apsidal-motion rates predicted 
by the classical formula deviate from the rates determined within the framework of the theory of dynamic tides 
due to the effects of the compressibility of the stellar fluid and du e to resonances of dynamic tides with free 
oscillation modes of the component stars jSmeyers fc Willems 20o"l ). As the star evolves on the main sequence, 
the deviations caused by the compressibility of the stellar fluid increase with increasing radius of the star. The 
additional deviations caused by the resonances are largest near the end of the core-hydrogen burning phase. Both 
of these deviations increase with increasing values of the orbital eccentricity. 
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1. Introduction 

In close binary systems of stars, each star experiences the 
time-dependent tidal force exerted by its companion. The 
distortion caused by the tidal force leads to a perturbation 
of the external gravitational field, which in its turn causes 
secular variations of the position of the periastron in the 
relative orbit of the companion. 

The most commonly used formula for the rate of secu- 
lar apsidal motion due to the tidal deforma tion of the com- 
ponents of a close bi nary was deri ved by [Cowling (1938 ) 
and subsequently by ^terne (1939 ) under the assumption 
that the orbital period is long in comparison with the free 
harmonic periods of the component stars. A generalisation 
of this classical formula taking into account the effects of 
an inclination of the equatorial planes of the stars with 
respect to the orbital plane of the binary was derived by 



centricities. When the orbital period is shorter, deviations 
arise due to the larger role of the stellar compressibility at 
higher forcing frequencies and due to the effects of reso- 
nances between dynamic tides and free oscillation modes. 
From the application of their results to 5 Mq, 10 Mq, and 



20 Mq zero-age main sequence stars, Smeyers & Willemi 



Kopal (1959| ). 

In a recent investigation, pmeyers" 



Willems (2001) 



(hereafter referred to as Paper I) reconsidered the deriva- 
tion of the classical apsidal-motion formula in the light of 
the theory of dynamic tides. The authors concluded that, 
for sufficiently long orbital periods, the tide raised by the 
companion can be approximated at each instant by an 
appropriate linear combination of three static tides. For 
binaries with sufficiently long orbital periods, the classical 
formula then turns out to be valid up to high orbital ec- 



found the relative differences between the classical formula 
and the formula established within the framework of the 
theory of dynamic tides to be inversely proportional to the 
square of the orbital period. The relative differences also 
turned out to be larger for a model with a larger mass. 

In this paper, we investigate the validity of the classical 
formula for the rate of secular apsidal motion in binaries 
with an evolved main-sequence star. Deviations between 
the classical formula and the formula established in the 
framework of the theory of dynamic tides may be expected 
to vary with the evolutionary stage of the star due to the 
changes in the central condensation and due to the growth 
of the stellar radius. 

The plan of the paper is as follows. In Sects. 2 and 3, 
we present the basic assumptions and the ingredients nec- 
essary for the determination of the rate of secular apsidal 
motion in binaries with eccentric orbits. In Sect. 4, we 
briefly recall the classical formula for the rate of secular 
apsidal motion and the corresponding formula established 
within the framework of the theory of dynamic tides. In 
Sect. 5, we examine the influence of stellar evolution on 
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the deviations between the apsidal-motion rates predicted 
by the classical formula and the rates given by the theory 
of dynamic tides for binaries with shorter orbital periods. 
Sect. 6 is devoted to concluding remarks. 



2. Basic assumptions 

Consider a close binary system of stars that are orbiting 
around each other in an unvarying Keplerian orbit with 
semi-major axis a and orbital eccentricity e. The first star, 
with mass Mi, is rotating uniformly around an axis per- 
pendicular to the orbital plane with an angular velocity 
which is assumed to be small in comparison to the inverse 
of the star's dynamical time scale: 



given by 



1 



Tdyn 



GMi 



1/2 



(1) 



Here G is the Newtonian constant of gravitation and i?i 
the radius of the undistorted spherically symmetric equi- 
librium star. We treat the companion star, with mass M2, 
as a point mass. 

Let Cixyz be an orthogonal right-handed frame of ref- 
erence whose origin coincides with the mass centre Ci of 
the tidally distorted star. The z-axis is normal to the or- 
bital plane and is oriented in the sense of the binary's or- 
bital angular momentum. The x- and y-axes are assumed 
to be corotating with the star. With respect to this frame 
of reference, we introduce a system of spherical coordi- 
nates r — (r, 0, (p) so that at each instant t the position of 
the companion is given by 

r = u, e = 7r/2, = w - f7 i, (2) 

where u is the radial distance and v the true anomaly of 
the companion in its relative orbit. 

The tidal force exerted by the companion is derived 
from the tide-generating potential eTW{r,t), where et 
is a small dimensionless parameter defined as 



Ri 



M2 
Ml' 



(3) 



The parameter corresponds to the ratio of the tidal force 
to the gravity at the star's equator. 



Following Polfliet & Smeyers (1990), we decompose the 
tide-generating potential in terms of unnormalised spher- 
ical harmonics y/"(0, (/)) and in Fourier series in terms of 
multiples of the companion's mean motion n = 27r/Porb, 
where Poib is the orbital period. The decomposition takes 
the form 

ErWir,t) = -Er^t E E 

e=2 m=-e fc=-oo 



Ri 



Ci,m,k \ I Ye' 



exp [i {aT t — krir) 



(4) 



where ax — k n+m is the forcing angular frequency with 
respect to the corotating frame of reference, r is a time of 
periastron passage, and the ce^rn,k are Fourier coefficients 



,|m| 



{£ + \m\y. 



(0) 



Ri 



1 

- {1 + e cosvY~'^ cos{kM + 
I" Jo 



(l-e2) 
m v) dv. 



e-1/2 



(5) 



In the latter expression, is an associated Legendre 

polynomial of the first kind, and M is the mean anomaly 
of the companion. 

The Fourier coefhcients ci^rn,k are equal to zero for odd 
values of ^-|- |to| and decrease with increasing values of the 
multiple k of the companion's mean motion. The decrease 
is slower for higher orbital eccentricities, so that the num- 
ber of terms that has to be taken into account in the ex- 
pansion of the tide-generating potential given by Eq. (jij) 
increases with increasing values of the orbital eccentricity. 

It follows that, for each degree £ of the spherical har- 
monics Y™{9, (j)), the tide-generating potential consists of 
a term giving rise to a static tide (ctt = 0), and an infi- 
nite number of terms giving rise to partial dynamic tides 
(ctt 7^ 0). The static tides are not to be assimilated with 
equilibrium tides which occur when the orbit of the com- 
panion is circular and the rotation of the star is synchro- 
nised with the orbital motion of the companion. 

3. Tlie equations governing forced isentropic 
oscillations of a spherically symmetric star 

When the effects of the Coriolis force and the centrifu- 
gal force are neglected, the uniformly rotating star can be 
considered to be spherically symmetric. In the isentropic 
approximation, a partial dynamic tide generated by a sin- 
gle term in the expansion of the tide-generating potential 
is then governed by the homogeneous fourth-order system 
of differential equations 



dr 

drjT _ 
dr 
1 d 

j,2 fj^j. 



9 



,d-^i 



dr 



£{£ + l) 



9 

^1) 



1 



Vt 



■qjr. 



Vt 



^1 



*T, (6) 



(7) 



= AnGp 



— 

9 



1 



— (erf, riT - ^'3 



(8) 



In these equations, Ct(^) and r]x{r) are the radial parts of 
the radial and the transverse component of the tidal dis- 
placement with respect to the local coordinate basis d/dr, 
d/dd, d/d(j), and '^ri^) is the sum of the radial parts of 
the tide-generating potential and the Eulerian perturba- 
tion of the star's gravitational potential. Furthermore, p is 
the density, g the gravity, the square of the isentropic 
sound velocity, and N'^ the square of the Brunt- Vaisala 
frequency. 

The solutions of the system of Eqs. - (||) must 
satisfy boundary conditions at the star's centre and at 
the star's surface. At r = 0, the radial component of the 
tidal displacement must remain finite. At r = Ri, the 
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Lagrangian perturbation of the pressure must vanish, and 
the continuity of the gravitational potential and its gradi- 
ent requires that 



dr 



i + 1 
Ri 



•-(2^+1) 



GM i 



(9) 



Because of the non-homogeneous term in the right-hand 
member of this boundary condition, the solutions of the 
system of Eqs. - (H) are proportional to the product 

ET Ce,m,k- 

In the case of static tides, the forcing angular frequency 
(Tt is equal to zero and the radial component of the tidal 
displacement is determined by the homogeneous second- 
order differential equation 



T;0 



2(1 ^ + 1 \ d^T;0 

J] dr r ) 



T:0 



0.(10) 



dr^ ' \g dr ' r J dr 
This equation is equivalent to the equation of Clairaut 
which is usually derived within the framework of the the- 
ory of equilibrium tides in the limiting case of an infinite 



orbital period ( Sterne 1939|) . 

The solution of Eq. ( |lO| ) must remain finite at ?' = 
and satisfy the boundary condition 



(eT;o)„ =eT (2^+1) 0,0- (11) 



V dr J rti 

For the remainder of the paper, we restrict the expan- 
sion of the tide-generating potential given by Eq. % to 
the terms associated with £ = 2, which are dominant. 



4. Apsidal motion due to the tidal deformation of 
a component of a close binary 

The rate of secular apsidal motion due to the tidal dis- 
tortion of the components of a close binary is usually de- 
termined by means of the classical formula established in 
the assumption that the orbital period is long in compar- 
ison to the periods of the free oscillation modes of the 
component stars (Cowling 1938, Sterne 1939). At each in- 
stant, the tide generated by the companion can then be 
approximated by a linear combination of three static tides 
(Paper I). The resulting formula takes the form 

5 



dTU 

~dt 



a 



Mo 2tt 



k2 15 / (e^ 



classical \ " / ^'^1 ^orb 

where w is the longitude of the periastron, and 



-5 



1 



1 



(12) 



(13) 



The constant fc2 is known as the apsidal-motion con- 
stant and depends on the mass concentration of the tidally 
distorted star. The constant is determined as 
'd h\{£_T-o/r) 



d Ini 



Ri 



d Inr 



(14) 



Because of the logarithmic derivatives in the numerator 
and the denominator of this definition, the apsidal-motion 
constant is independent of the product Et C2,o,o appearing 
in the boundary condition given by Eq. (pT|). 

Within the framework of the theory of dynamic tides, 
the rate of secular apsidal motion is determined by adding 
the contributions to the apsidal motion stemming from the 
various partial static and dynamic tides which are gener- 
ated by the individual terms in the expansion of the tide- 
generating potential given by Eq. (|^). The resulting rate 
of secular change of the longitude of the periastron is given 

by 



dTu\ 



Ri 

a 



Ah 2 7r 



2 fc2 G2.o,o + 4 



Ml Porb 



E 

k=l 



(15) 



where the G2.m,k are functions of the orbital eccentricity 
defined as 

G2,m,k{e) = — ^— -2C2,m,fc^'i'"'(0) 



e {l-e^Y 



3 / (1 -I- e cos w)^ cos(m w -I- fc M) cos w 





-m I (1 + ecosw) (2 -f ecosu) 
/o 



sin(m V + k M) sin v dv 



(16) 



( ISmeyers ct al. 1998| ). 

The coefficients F2^m,k render the star's responses to 
the various forcing angular frequencies ctt and are related 
to the total perturbation of the gravitational potential at 
the star's surface as 

\_G Ml £t C2,m,k 

In analogy with the apsidal-motion constant ^2, the co- 
efficients F2,m,k are independent of the product £t C2,m k 
which appears in the boundary condition given by Eq. (g) . 

In the limiting case of long orbital and rotational pe- 
riods, the classical formula for the rate of secular apsi- 
dal motion agrees with the formula established within the 
framework of the theory of dynamic tides up to high or- 
bital eccentricities ( Smeyers et al. 199§| ). The agreement 
rests on the property that, in this limiting case, all forcing 
angular frequencies (tt tend to zero, so that the coeffi- 
cients F2.m,k, with m = —2, 0, 2 and k = 1,2,3,..., tend 
to the apsidal-motion constant k2. 

For binaries with shorter orbital periods, deviations 
arise due to the increasing role of the stellar compress- 
ibility at higher forcing frequencies and due to resonances 
of dynamic tides with lower-order (/"'"-modes of the com- 
ponent stars (Paper I). The extent of the deviations can 
be expected to depend on the evolutionary stage of the 
star due to the increase of the stellar radius and the re- 
distribution of the mass as the star evolves on the main 
sequence. 
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Fig. 1. Evolutionary track of the 5 Mq star with initial 
chemical composition X = 0.70 and Z = 0.02. The labels 
correspond to the numbers of the models considered in 
our investigation. 



5. Deviations of the classical apsidal-motion 
formula for binaries with an evolved 
main-sequence star 

5.1. The models 

We determined the deviations of the classical formula for 
the rate of secular apsidal motion in binaries with shorter 
orbital periods for a sequence of models of a 5 Mq star 
with initial chemical composition X = 0.70 and Z = 0.02. 
The evolutionary track of the star in the Hertzsprung- 
Russell diagram is displayed in Fig. |l|. The labels in the 
figure correspond to the numbers of the models considered 
in our investigation, which range from the last stages of 
the C^^ N^** reaction (model 20) to the phase where 
hydrogen is exhausted in the core (model 96). 

The variation of the central hydrogen abundance Xc, 
the apsidal-motion constant k2, and the radius Ri are 
displayed in Fig. || as a function of the age of the star. 
The central hydrogen abundance Xc decreases as the star 
evolves along the main sequence due to the conversion of 
hydrogen into helium via the CNO cycle. In addition, the 
conversion of hydrogen into helium increases the central 
condensation of the star, so that the value of the apsidal- 
motion constant k2 decreases. At the end of the main se- 
quence, when the star has reached an age of 8.4 x 10'' 
years, the central hydrogen abundance has dropped to 
0.001 and the apsidal-motion constant has decreased by 
a factor of two. The radius of the star expands during 
the major part of the evolution on the main sequence. It 
reaches a maximum value for model 72 and subsequently 



Fig. 2. The variations of the central hydrogen abundance 
Xc, the apsidal-motion constant A;2, and the radius Ri as 
a function of the age of the 5 Mq star. 

Table 1. Properties of the models of the SM© star con- 
sidered in this investigation. 



Model Nr. 


age (yrs) 


Ri (km) 




k2 


20 


3.03 10^ 


1.93 10** 


0.699 


0.00741 


52 


5.65 10^ 


2.66 10** 


0.358 


0.00505 


72 


8.13 10^ 


3.63 10^ 


0.060 


0.00369 


96 


8.44 10^ 


3.29 10** 


0.001 


0.00348 



decreases again during the final stages of core-hydrogen 
burning. The ages, the radii, the central hydrogen abun- 
dances, and the apsidal-motion constants of the models 
used in our investigation are listed in Table ^ 

Due to the growth of the radius with the age of the 
star, care needs to be taken that only orbital periods are 
considered for which the radius of the star remains smaller 
than the radius of its Roche lobe. For binaries with circular 
orbits, the star's Roche-lobe radius Rl.i can be approxi- 



mated by means of Eggleton's (1983) fitting formula 
i?L4 0.49g-2/3 



(18) 



a 0.6g-2/3 + ln(l + g-i/3)' 

where q = M2/M1. The relative error of the formula is 
smaller than 2% for all values of the mass ratio in the 
range oi < q < 00. 

The variation of the Roche-lobe radius of a 5 Mq star, 
expressed in solar radii, is presented in Fig. ^ as a function 
of the orbital period for the mass ratios q = 0.2, 0.5, 1.0. 
The solid lines in the figure correspond to the Roche-lobe 
radii for binaries with a circular orbit. The dotted hori- 
zontal lines represent the stellar radii at the evolutionary 
stages represented by the models 20, 52, 72, and 96. 
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Fig. 3. Variation of the Roche-lobe radius of a 5 M0 star 
in the case of the orbital eccentricities e = 0.0,0.25,0.5 
and the mass ratios q — 0.2, 0.5, 1.0. The dotted horizontal 
lines represent the radius of the star for models 20, 52, 72, 
and 96. 



In the case of a binary with an eccentric orbit, the 
Roche-lobe radius can be expected to be smallest when 
the companion is located in the periastron of its relative 
orbit. Therefore, we also determined the variations of the 
Roche-lobe radius for binaries with an eccentric orbit by 
substituting the periastron distance Tp = a(l — e) for the 
semi-major axis a in Eq. (|l8|). The resulting Roche- lobe 
radii are represented by the long- and the short-dashed 
lines in Fig. || for the orbital eccentricities e = 0.25 and 
e = 0.5, respectively. 

The Roche-lobe radii increase monotonically with in- 
creasing values of the orbital period and are larger for 
smaller values of the mass ratio. The restrictions on the 
orbital period imposed by the requirement that the star fit 
within its Roche-lobe are most stringent for higher orbital 
eccentricities. In the case of model 72, e.g., the orbital 
period must be larger than 3 days when the eccentricity 
takes the value e = 0.25 and larger than 5.5 days when 
the eccentricity takes the value e = 0.5. 

5.2. Deviations of the classical formula 



An ap]3ropriate quantity to examine the deviations be- 
tween the rate of secular apsidal motion given by the 
classical formula and the corresponding rate determined 
within the framework of the theory of dynamic tides is the 
relative difference 



Table 2. The proportionality factor A for models 20, 52, 
72, and 96 of the 5 Mq star and the orbital eccentricities 
e = 0.25 and e = 0.5. 



Model Nr. 


Ae=0.25 


Ae=0.5 


20 


0.13 


0.52 


52 


0.27 


1.1 


72 


0.88 


3.2 


96 


0.68 


2.3 



This relative difference is independent of the product 
(i?i/a)'(M2/Mi)(2 7r/Porb), 

so that it depends on the orbital period only through the 
coefficients F2^m,k- 

We determined the relative differences A as a function 
of the orbital period for the evolutionary models 20, 52, 
72, and 96 of the ^Mq star, and the orbital eccentricities 
e = 0.25 and e = 0.5. The shortest orbital periods consid- 
ered for each model are chosen such that the radius of the 
star is sufficiently small in comparison to the radius of its 
Roche lobe (Sect. 



5.1 



For the determination of the forc- 
ing angular frequencies ctt with respect to the corotating 
frame of reference, we adopted the low angular velocity of 
rotation O = 0.01 n, so that the companion's orbital pe- 
riod with respect to the corotating frame remains short. 
In the case of a resonance of a dynamic tide with a free os- 
cillation mode, we limit the determination of A to orbital 
periods for which the relative difference between the forc- 
ing angular frequency of the tide and the eigenfrequency 
of the oscillation mode is larger than 0.1 £t (see Paper I). 
The resulting variations of the relative difference A are 
represented by the solid lines in the left- and right-hand 
panel of Fig. ^. 

The relative differences A are mostly negative so that 
for binaries with shorter orbital periods the classical for- 
mula yields somewhat too small values for the rate of sec- 
ular apsidal motion, and thus somewhat too long apsidal- 
motion periods. The peaks observed at shorter orbital pe- 
riods are caused by resonances of dynamic tides with free 
oscillation modes (7+ of the tidally distorted star. These 
peaks are superposed on a basic curve which represents 
the systematic deviations caused by the increasing role of 
the stellar compressibility at shorter orbital periods. 

The systematic relative differences A in the absence of 
resonances can be approximated by a formula of the form 
Ataxic = -A r-2, (20) 

where A is a constant. For the determination of A we used 



a linear least-squares procedure described by Press et al 



A = 



{dw/dt)^^^ 



(19) 



(1992) and fitted Eq. ( |20| ) to the relative differences A 
displayed in Fig. ^ The resulting values of A are listed 
in Table || in the supposition that the orbital period is 
expressed in days. The associated basic curves determined 
by means of Eq. ( p0| ) are represented by the dashed lines 
in Fig. ^. 

The downwards slope of the basic curves is larger for a 
model representing a more evolved stage of the 5 Mq star. 
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Fig. 4. The relative differences A for models 20, 52, 72, and 96 of the 5 Mq star in the case of the orbital eccentricities 
e = 0.25 (left-hand panel) and e = 0.5 (right-hand panel). The solid lines represent the relative differences in the 
presence of resonances of dynamic tides with free oscillation modes g~^. The dashed lines correspond to the relative 
differences in the absence of resonances. 



The slope is largest for model 72, which is the model with 
the largest radius, and decreases again during the final 
stages of core-hydrogen burning. The underlying reason 
can be identified by expressing the forcing angular fre- 
quencies <tt in units of the inverse of the star's dynamical 
time scale: 



(It 



\ 1/2 



these resonances the relative differences between the rates 
of secular apsidal motion predicted by the classical for- 
mula and the rates given by the theory of dynamic tides 
can amount to as much as a few tens of percents, espe- 
cially for higher orbital eccentricities. This conclusion is 
in agreement with the conclusions reached by Papaloizou 
fc Pringle (198C| ) and [Quataert et al. (19961) . 



(21) 



For a model with a constant mass, it then follows that 
the dimensionless forcing angular frequencies increase 
with increasing values of the stellar radius. The systematic 
deviations caused by the compressibility of the stellar fluid 
will therefore be larger for a model with a larger radius. In 
the case of the orbital eccentricity e = 0.5 and the orbital 
period of 6 days, e.g., the relative differences A in the 
absence of resonances increase from 1.4% for model 20 to 
8.9% for model 72. 

In the case of a resonance of a dynamic tide with a 
free oscillation mode, the relative differences A may de- 
viate substantially from the systematic trend represented 
by the basic curves. The largest deviations occur for res- 
onances between the lower-order harmonics in the expan- 
sion of the tide-generating potential given by Eq. (||) and 
the lower-order g^-modes of the tidally distorted star. For 



For a given range of orbital periods and a fixed or- 
bital eccentricity, the deviations caused by the resonances 
increase as the star evolves on the main sequence. This be- 
haviour is related to the growing radiative envelope which 
causes the eigenfrequencies of a more evolved stellar model 
to be smaller and closer to each other than those of a less 
evolved stellar model, so that resonances with lower-order 
(7 "'"-modes become increasingly important with increasing 
age of the star. A similar increase in the effectiveness of 
tidal forcing during the late stages of core hydrogen burn- 
ing was found by Savonije fc Papaloizou (1984 ) in an in- 
vestigation on the spin-up of neutron stars and the circu- 
larization of orbits in massive X-ray binaries. 



6. Concluding remarks 

In this paper, we have examined the influence of stellar 
evolution on the deviations between the rates of secular 
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apsidal motion predicted by the classical formula due to 
Cowling (193^ ) and Sterne (1939| ), and the rates of sec- 



ular apsidal motion determined in the framework of the 
theory of dynamic tides. The deviations arise for binaries 
with shorter orbital periods due to the growing role of the 
stellar compressibility at higher forcing frequencies and 
due to resonances of dynamic tides with the free oscilla- 



tion modes of the component stars (Smeyers & Willem; 
200l|)~| 

The extent of the deviations caused by the compress- 
ibility of the stellar fluid depends on the evolutionary stage 
of the star primarily through the evolution of the radius 
and the associated change of the star's dynamical time 
scale. As the star evolves on the main sequence, the ra- 
dius and the dynamical time scale increase so that the 
orbital period effectively becomes shorter in comparison 
to the star's dynamical time scale. The deviations due to 
the compressibility of the stellar fluid are therefore larger 
for a model with a larger radius. 

Besides varying the systematic deviations caused by 
the compressibility of the stellar fluid, the evolution of the 
star also affects the deviations caused by the resonances 
of dynamic tides with free oscillation modes of the compo- 
nent stars. The latter deviations are superposed upon the 
deviations caused by the stellar compressibility and can 
amount to as much as ten percent and more, especially 
for higher orbital eccentricities. In the case of the 5 Mq 
star considered, the effects of the resonances are larger for 
a model near the end of core-hydrogen burning than for a 
model near the zero-age main sequence. 

Acknowledgements. Bart Willems acknowledges the support of 
PPARC grant PPA/G/S/1999/00127. 
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